Spontaneous breaking of four-fold rotational symmetry in two-dimensional electronic 
systems explained as a continuous topological transition 
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The Fermi liquid approach is applied to the problem of spontaneous violation of the four-fold 
rotational point-group symmetry (C4) in strongly correlated two-dimensional electronic systems on 
a square lattice. The symmetry breaking is traced to the existence of a topological phase transition. 
This continuous transition is triggered when the Fermi line, driven by the quasiparticle interactions, 
reaches the van Hove saddle points, where the group velocity vanishes and the density of states 
becomes singular. An unconventional Fermi liquid emerges beyond the implicated quantum critical 
point. 
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The breaking of fundamental symmetries in ground 
states of strongly correlated two-dimensional (2D) elec- 
tron systems^"— remains one of the most intensely 
debated topics in low-temperature condensed matter 
physics. Kivelson, Fradkin, and Emery 7 were the first 
to discuss the case of nematic phase transitions, well 
before relevant experimental data was obtained. Some- 
what later, Yamase and Kohnc 8 (within t — J model) and 
Halboth and Metzner 9 (within the Hubbard model) at- 
tributed the breaking of four-fold symmetry to violation 
of a Pomeranchuk stability condition 10 associated with 
antiferromagnetic fluctuations. 

Subsequently, much theoretical work has been aimed 
at elucidating salient features of this phenomenon, pri- 
marily within mean-field theoryii - — It is instructive 
to recognize that the approach taken in these ef- 
forts bears a striking resemblance to that employed by 
Belyaev fifty years ago to describe quadrupole defor- 
mation in atomic nuclei. 18 To determine the critical 
point at which the spherical shape becomes unstable 
and calculate the nuclear deformation beyond this point, 
he introduced an effective Hamiltonian with separable 
quadrupole-quadrupole Q1Q2 interaction. Analogously, 
for two-dimensional tetragonal electronic systems, a sep- 
arable interaction <^2(pi)^2(P2) with order parameter 
d2(p x ,Py) = cos p x — cos p y is adopted in the mean- field 
treatments of the breakdown of C4 symmetry, momenta 
being measured in units of the inverse lattice constant. 
However, such an effective Hamiltonian with separable 
interaction is appropriate only in the channel where sym- 
metry breaking occurs. Moreover, even in this channel a 
mean- field approach may be inadequate, as exhibited for 
example in the prediction of a first-order phase transition 
in the case of violation of point-group symmetries on a 
square lattice. 19 

Burdened with variety of inconsistencies, the mean- 
field description of nuclear deformation was superseded 
many years ago by the more sophisticated Fermi-liquid 
(FL) approach^ Following this successful precedent, we 
work within the FL framework to obtain a better under- 



standing of C^-symmetry breaking in electron systems 
on a 2D square lattice. Intensive numerical calculations 
assuming a finite-range exchange interaction, supported 
by complementary analysis of a simplified model, disclose 
unexpected features of the phenomenon. In contrast to 
the description given by mean-field theory, we find that 
the breakdown of C4 symmetry is associated with a topo- 
logical phase transition that occurs under conditions that 
allow the Fermi line, calculated within FL theory, to reach 
the van Hove saddle points (0, 7r), (7r, 0), (0, — 7r), (— 7r, 0). 

Consideration of topological transitions dates back to 
an article by I. M. Lifshitz 21 , in which the form of 
the single-particle spectrum e(p) was assumed to be 
known. However, within FL theory e(p) is itself a 
functional of the quasiparticle momentum distribution 
n(p) = [1 + exp ((e(p) — /i)/T)] _1 . Accordingly, self- 
consistent inclusion of the interactions between quasi- 
particles may lead to unforeseen types of the topological 
transitions^"— We find just such a case in the problem 
of (^-symmetry violation. 

Stated simply, topological transitions in correlated 
Fermi systems are signaled (at zero temperature) by a 
change of the number of roots of equation 



e(p,n F ) = /i, 



(1) 



where rip is the Fermi step distribution and \i is the 
chemical potential. For a thorough development of the 
concept, see the review by Volovik^ - Throughout, we ad- 
here to his rigorous quantitative definition of topological 
phase transitions, as distinguished from looser notions 
such as transitions between large and small Fermi sur- 
faces that are also prevalent in the literature. 

Analysis of topological phase transitions in fermionic 
systems is greatly facilitated by the absence of criti- 
cal fluctuations of any order parameter at the transi- 
tion point and its vicinity, meaning that the Landau- 
Migdal quasiparticle picture retains its validity. Thus, 
the physical many-fermion system may be viewed as a 
system of interacting quasiparticles, and C4 symmetry 
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violation can be investigated using the fundamental FL 
relation^^ 

<9e(p) <9e (p) I f dn(pi) 

-aT = ^^ + 2 Tr y ^^(P'Pi)-^-^i- ( 2 ) 

In this relation, dv = dp x dp y / (2tt) 2 is the volume element 
of 2D momentum space, eo(p) is the bare single-particle 
spectrum, and J 7 is a phenomenological interaction func- 
tion depending only on the momenta p and pi of the 
colliding quasiparticles. 

Our goal is to analyze the impact of antiferromagnetic 
fluctuations on electron spectra calculated using Eq. (j2j). 
Taking account of these fluctuations in the interaction 
function T presents little difficulty in the regime far from 
the antiferromagnetic phase transition, since the fluctua- 
tion exchange can be treated within the Ornstein-Zernike 
approximation. The part of T responsible for the ex- 
change is then given by 

K^siP, Pi) = AV^ov [(p-Pi-Q) 2 + r 2 ] _1 , (3) 

where the constant A represents the spin-fluctuation ver- 
tex and Q = (7r,7r) the antiferromagnetic wave vector, 
while £ denotes the correlation radius. Result ([3]) relies 
on the fact that the interaction function T coincides with 
a specific static limit of the quasiparticle scattering am- 
plitude whose initial and final energies are on the Fermi 
surface, such that this quantity is energy- and frequency- 
independent^^ 

Inserting Eq. (|3]) into Eq. (|2]) and evaluating the spin- 
fluctuation contribution aided by the identity 2cr a f3cr 1 5 — 
3^(5^7/3 — GolSGiP, one arrives at 

e(p)=e (p) + /a / (p _ p ;^ ) )2+c _ 2 ^ 1 , (4) 

where f a = 3A 2 /2. The chemical potential /i, being con- 
stant along the Fermi line, is determined by the normal- 
ization condition j n(p) 2dv = p, where factor 2 assumes 
summation over two spin projections. 

Numerical solution of the 2D nonlinear integral equa- 
tion (jlj) is extremely time-consuming since it is necessary 
to compute to high precision to rule out spurious signals 
of broken symmetry. 

Calculations have been performed in the case of an 
open Fermi surface, assuming a tight-binding spectrum 

eo(p) = -2 to (cosp^+cosp^+^i cos p x cos py-/i, (5) 

with the ratio ti/to of input parameters to and t\ taken to 
ensure a rather small distance between the tight-binding 
Fermi line and the saddle points. A finite-range interac- 
tion function 

/(q)=/a[(q-Q) 2 +r 2 ] _1 (6) 

is adopted, with £ = 30. Salient results are reported 
for strategically chosen values f a = 1.0 and 1.5 of the 
coupling parameter (in units of 2to), and at temperatures 
T = 10~ 2 and 10~ 4 (also in units of 2t ). 



The numerical calculations, as represented here in 
Figs. [l]-[3j reveal some remarkable features of the quasi- 
particle rearrangement responsible for violation of C4 
symmetry. One conspicuous feature, well documented 
in the figures, is that only those quasiparticles residing 
in domains close to the saddle points are noticeably af- 
fected by the inclusion of antiferromagnetic fluctuations 
in the FL treatment. 

To understand of the onset of symmetry breaking, it 
is most instructive to track the distance between neigh- 
boring points where the Fermi line crosses the border of 
the Brillouin zone. This distance is found to shrink as 
the coupling constant f a is increased toward a critical 
value lying between 1.0 and 1.5. Breakdown of C4 sym- 
metry presumably occurs for a critical coupling f ac at 
which the two crossing points merge with one another as 
they embrace and converge upon the nearby saddle point. 
As seen in Fig. [TJ the Fermi line (in red) calculated at 
temperature T = 10 -4 (effectively zero) for a coupling 
constant f a = 1.5 lying beyond the critical point does 
indeed violate x — y symmetry; the corresponding Fermi 
line for f a = 1.0 (in green) does not. Another feature 
worthy of note is the effect of temperature in suppressing 
the symmetry breaking phenomenon. Comparison of the 
deviations of the Fermi line calculated at f a = 1.5 and 
T = 10 -4 (in red) from the symmetry-preserving Fermi 
line obtained at f a = 1.5 and T = 10 -2 (in blue) demon- 
strates that C4 symmetry can be restored by elevating 
the temperature. 

Results from calculations of the magnitude v(p) = 
\de(p)/dp\ of the group velocity along the Fermi line are 
plotted in Fig. [2j These results demonstrate that the 
impact of antiferromagnetic correlations, as described by 
Eq. (j3]), is only significant for quasiparticles in momen- 
tum domains adjacent to the saddle points. The small 
gap between the values for v given by the bare tight- 
binding model and by the Fermi-liquid-theory treatment, 
seen in domains away from the saddle points, is due pri- 
marily to a shifting of the location of the Fermi line trig- 
gered by the antiferromagnetic correlations. On the other 
hand, the group velocity v evaluated at f a = 1.5 and 
T = 10 -4 (red line) is seen to acquire an x — ^-anisotropic 
component close to the saddle points. 

Fig. [3] presents results from calculations of electron 
spectra at T = 10 -4 in the direction perpendicular to 
the Fermi line. Far from the saddle points, the impact of 
antiferromagnetic fluctuations is again found to be mi- 
nor, but in their vicinity the effects are very strong. In 
particular, the particle and hole spectra cease to be alike; 
the average slope of the hole spectrum noticeably exceeds 
that of the particle spectrum. One might attribute this 
difference to the variation of Fermi-line contributions to 
Eq. ([2|) associated with a turning point emerging in the 
trajectory of the Fermi line at the critical point. From all 
the results discussed above, we infer that the description 
of the rearrangement of the ground state in terms of a 
single di parameter is a poor approximation. 

To further analyze and interpret the results obtained 
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FIG. 1: Fermi lines computed for the Fermi- liquid-theory 
model assuming bare tight-binding spectrum © with ti/to = 
0.45 and finite range interaction (|6)) with £ = 30. Panel (a): 
Results for f a = 1.5 and T = 10 -4 (both in units of 2 to). 
Thick solid line (in red): one of two identical solutions with 
spontaneously broken Ca symmetry. Only the first quadrant 
of the Brillouin zone is drawn, since neither p x — > —p x nor 
p y — >• — p y reflection symmetry is broken. Thin solid lines (in 
black) : Fermi lines for the bare tight-binding spectrum e° and 
its counterpart. Panels (b) and (c): Two shaded squares ad- 
jacent to the saddle points (0,7r) and (7r,0) present in panel 
(a) are magnified. The Fermi- line solution with broken Ca 
symmetry appearing in panel (a) (red line) is drawn together 
with two x — ^-symmetrical solutions corresponding respec- 
tively to f a = 1.0 and T = 10 -4 (green line), and to f a = 1.5 
and T= 10" 2 (blue line). 



numerically for the finite-range interaction (j6j) and bare 
tight-binding spectrum ([5]), we simplify the interaction 
function in the manner of Ref. [27|, replacing the exchange 
term (j3j by an infinite-range form 

/(q) = (27r) 2 /o<5(q-Q), (7) 

with coupling constant /q. Eq. (j4| is then replaced by 27 

e(p, T) = e (p) + /on(e(p + Q, T)). (8) 

This treatment is analogous to that adopted by 
Nozieres 24 in a study of non-FL behavior of strongly cor- 
related Fermi systems for the case where forward scatter- 
ing is dominant. Eq. (j8j) can be derived within a standard 
variational procedure based on the formula 27 



E 



I 
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Cp"(P)+ o/0 n (PMP + Q) 



2dv 



(9) 




FIG. 2: Group-velocity magnitudes v — \de(p)/dp\ (in units 
of 2to), evaluated along the Fermi line as a function of the an- 
gle (p defined in the inset, for different single-particle spectra 
e(p). Results are shown for the bare tight-binding model with 
the same parameter choice as in Fig. [T] (brown line) and for 
the Fermi- liquid-theory model of Fig. [1] at f a = 1.0, T = 10 -4 
(green line); f a = 1.5, T = 10 -4 (red line); and f a = 1.5, 
T — 10 -2 (blue line). Broken C4 symmetry of the red curve 
with respect to x — y exchange is manifested by its different 
behavior in the two shaded areas close to the saddle points. 



Eq. flSJ is conveniently rewritten as a system of two equa- 
tions 



e i = e i 

^2 



■fon(e 2 ), 
4 + /oro(ei), 



(10) 



for the energy E of the model quasiparticle system. 



where e\ = e(pi) — ji and 62 = e(pi + Q) — /i, while 
e? = e (p) and = e (p + Q). 

In the earlier work, 27 a graphical procedure was intro- 
duced to solve the set ([TO]) at T = 0. Three different 
solutions were found. One of these corresponds to an ex- 
ceptional, non-FL stated 2 - exhibiting a flat single-particle 
spectrum. In the absence of pairing correlations, this so- 
lution turns out to be disfavored energetically relative to 
the other two solutions, which possess identical FL-like 
properties. 

Focusing on the properties of the latter two solutions, 
we observe that at T — the associated rearrangement 
of the initial Landau state can occur only in those 2D 
systems in which hot spots 36 exist — points situated on 
the Fermi line and connected by the vector Q. In fact, 
for systems with small quasiparticle filling, the product 
n(p)n(p + Q) vanishes for any momentum p; hence the 
ground-state energy is independent of the coupling con- 
stant fo. The same is true in the case of small quasihole 
filling. 

In systems having hot spots, the rearrangement occurs 
due to breaking of quasiparticle pairs occupying single- 
particle states with momenta p and p + Q. The cor- 
responding domain 1Z (the "reservoir") consists of four 
quasi-rect angles, each adjacent to one of the van Hove 



4 




71 



FIG. 3: Single-particle spectra e(p) (in units of 2to) evaluated 
along two lines in the momentum plane (as indicated with 
corresponding color coding in the insets). The Fermi- liquid 
theory model used for Fig.[T]is applied at T=10 -4 with / a = 1.0 
(upper panel) and / a =1.5 (lower panel). Spectral curves are 
plotted versus the normal component p n of the momentum, 
measured with respect to the Fermi line. In the insets, the 
blue line coincides with the diagonal of the zone, while the 
red line crosses the relevant hot line. 



saddle points. Each of the four elements of 1Z is confined 
between the border of the Brillouin zone, the counter- 
part of the initial Fermi line, defined by the equation 
eo(p + Q) = and two segments of the Fermi line em- 
bracing the given saddle point. 

In the rearrangement being considered, the quasipar- 
ticles move out the domain 1Z to resettle in a region 
C where all pairs of single-particle states connected by 
the vector Q are empty. The region C comprises four 
"lenses," situated between neighboring hot spots and 
bounded by the initial Fermi line and its counterpart (see 
panel (a) of Fig. 2]). The transfer of one quasiparticle 
from TZ to C produces a gain in energy which is just the 
coupling constant fo minus the loss r of kinetic energy. 
The minimum loss r m i n occurs when a quasiparticle, va- 
cating a state in 1Z with momentum p, occupies in £ a 
state of lowest energy, given by the chemical potential, 
so that T m i n = fi — eo(p). Therefore the rearrangement is 
favorable provided eo(p) — M + /o > 0. In the resettlement 
process, the chemical potential /i, which coincides with 
the maximum quasiparticle energy in occupied states (in 
particular, in the lens region), evidently increases relative 
to its initial value fa. The quasiparticles that resettle to 
the lens region then possess almost the same effective 
mass as the noninteracting electrons. This conclusion is 
confirmed by the numerical calculations represented in 
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FIG. 4: Panel (a): Fermi line (black) and its counterpart 
(blue) for the bare tight-binding spectrum of Eq. ([5]) with 
ti/to = 0.45. The "reservoirs" 1Z (see text) are colored in 
green, and the lenses C, in light gray. The hot spots con- 
nected with each other by the vector Q are symbolized by 
red dots. Panel (b): Fermi line for the simplified Fermi- liquid- 
theory model based on the infinite-range interaction function 
© with / =0.4 (in units of 2to). Hot lines are drawn as red 
line. Fermi lines for the bare tight-binding spectrum eo(p) 
and for the same spectrum shifted by —fo are shown as dashed 
lines while their counterparts are drawn in blue. 



Fig. El 

An alternative process involves transfer of the quasi- 
particle counterpart, which has momentum p+Q. In this 
case, the rearrangement occurs provided eo(p + Q) — fi + 
/ > 0. The choice between the two options is decided by 
comparison of the corresponding energies. The boundary 
at which one behavior gives way to the other is defined 
by the relation eo(p) = eo(p + Q). Since the straight line 
so defined is part of the new Fermi line, we infer that the 
rearrangement has converted the original, isolated hot 
spot into a continuous line of hot spots (see panel (b) of 

Fig. ED. 

The results obtained imply that quasiparticles are 
swept from a certain subdomain S of 1Z consisting of 
eight approximately trapezoidal strips. The boundaries 
of a given strip are traced on three sides by (respectively) 
the initial Fermi line, the border of the Brillouin zone, 
and a line geometrically similar to the initial Fermi line 
but shifted into the domain 1Z (see Fig. [4j). The strip's 
fourth side (red line) is just the hot line. The solution 
derived is self-consistent: any single-particle state with 
momentum p G S has its counterpart, with momentum 
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p + Q, located outside <S, and this state is occupied, so 
that Eq. (|5J) is fulfilled. Transparently, in this non-critical 
scenario, the new momentum distribution does not vio- 
late C4 symmetry. 

In the situation where C4 symmetry is violated in the 
rearrangement, the symmetry breaking occurs for a criti- 
cal value f c of /o , at which two segments of the Fermi line 
crossing the same boundary of the Brillouin zone merge 
at the saddle point. When this happens, the number of 
solutions of Eq. (pp) certainly drops, signaling a topolog- 
ical phase transition which, as readily seen, entails the 
breakdown of C4 symmetry. 

Suppose on the contrary that C4 symmetry is pre- 
served at /o > f c - Then all the saddle points must be 
emptied simultaneously, implying that every rearranged 
saddle point energy e s exceeds the chemical potential \i. 
But according to Eq. ((8]), the interaction contribution to 
e s vanishes when all the saddle points are emptied. Con- 
sequently, at /o > / c , the saddle-point energy e s must 
coincide with the corresponding bare value e^, which is 
lower than the initial chemical potential \±i. Thus, a con- 
tradiction is encountered. 

We are driven to the conclusion that the critical sit- 
uation giving rise to violation of C4 symmetry is one in 
which the Fermi line, calculated within FL theory, attains 
a saddle point. Since both components of the quasiparti- 
cle group velocity v(p s ) vanish at this critical point, the 
corresponding density of states must acquire a singular- 
ity, which implies that we are dealing with a quantum 
critical point (QCP). 

The contradiction is resolved beyond the QCP if only 
one of two neighboring saddle points is emptied, with the 
second remaining occupied — thereby breaking C4 sym- 
metry. As a point where the Fermi line crosses the p x 
axis moves away from the affected saddle point, its coun- 
terpart, shifted by the vector Q, slides along the border of 



the Brillouin zone, determining the boundary of the new 
filling. These conclusions drawn from analysis of the sim- 
ple infinite-range model are in agreement the findings of 
the numerical calculations based on the more elaborate 
model based on Eqs. (|4])-(j6]). 

To summarize: in addressing the problem of C4- 
symmetry violation, we have taken account of antiferro- 
magnetic fluctuations within a self-consistent Fermi liq- 
uid approach, employing an interaction function that is 
more realistic than the separable approximation assumed 
in mean-field treatments. We have demonstrated that in- 
clusion of the exchange interaction drives the calculated 
single-particle spectrum so as to shrink the distance be- 
tween saddle points and the Fermi line. When merging 
occurs, the electron group velocity vanishes at the points 
of mergence, because these points coincide with the sad- 
dle points. A quantum critical point (QCP) of a new type 
is thereby revealed, at which a topological phase transi- 
tion triggers the violation of C4 symmetry. Significantly, 
the transition is found to be continuous, in contrast to the 
first-order phase transition obtained in mean-field theory, 
where the corresponding QCP does not exist. Beyond the 
transition point, the group velocity becomes finite again. 
Thus, on one side of the QCP, the system behaves as 
conventional Landau Fermi liquid. On the other side, 
the electron liquid becomes an unconventional Fermi liq- 
uid because of the loss of four-fold symmetry. 
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